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Abstract 

We show that any pair of real symmetric tensors e and /j, can be realized as the 
effective electric permittivity and effective magnetic permeability of a metamaterial 
at a given fixed frequency. The construction starts with two extremely low loss 
metamaterials, with arbitrarily small microstructure, whose existence is ensured 
by the work of Bouchitte and Bourel and Bouchitte and Schweizer, one having at 
the given frequency a permittivity tensor with exactly one negative eigenvalue, and 
a positive permeability tensor, and the other having a positive permittivity tensor, 
and a permeability tensor having exactly one negative eigenvalue. To achieve the 
desired effective properties these materials are laminated together in a hierarchical 
multiple rank laminate structure, with widely separated length scales, and varying 
directions of lamination, but with the largest length scale still much shorter than 
the wavelengths and attenuation lengths in the macroscopic effective medium. 
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Introduction 

Isotropic materials with negative values of permittivity e and/or permeability fi (with 
low loss) have some fascinating properties. The colors of some stained glass windows are 
due to a resonance effect caused by metal spheres, which have a negative permittivity, 
embedded in the glass pp. Veselago [2] found that a slab of material with e — \i — — 1 at 
a given frequency would have a negative refractive index and act as a lens. Nicorovici, 
McPhedran and Milton [3] found that a cylindrical shell with permeability e — — 1 would, 
in the quasistatic limit and for TM waves, magnify the core material and produce a perfect 
image of a line dipole in the limit as the loss went to zero. Pendry [3] made the remarkable 
assertion that the Veselago lens would behave as a superlens, not limited by diffraction, 
providing perfect images of point sources, not just in the quasistatic limit, and now the 
validity of this assertion is beyond doubt provided one adds an infinitesimal loss to the 
lens: see, for example, [51 EJ [7J El [9] • Polarizable dipoles can be cloaked by such lenses [TUl 
EH H21 [131 E] ) an d larger objects can be cloaked by embedding the matching "antiobject" 
in the lens [15]. However any reasonable losses in the material can almost destroy these 
effects. From a practical viewpoint losses are always present and even extremely tiny 
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losses can prevent the functioning of the Pendry superlens if the wavelength is short 
compared with the lens thickness [SI [16] . However the negative refractive index property 
is robust: materials with a negative refractive index were manufactured by Shelby, Smith 
and Schultz [T7] . 

Other remarkable applications are associated with anisotropic materials having per- 
mittivity tensors e and permeability tensors fj, with e = /j,. It has long been known that 
the time harmonic Maxwell's equations are invariant under curvilinear coordinate trans- 
formations [18] and Dolin [19] realized that the transformed equations can be viewed 
as a solution of Maxwell's equations in a new material. Using such transformations 
empty space with e = // = 1 is transformed in a material with e(x) = //(x): thus Dolin 
recognized that certain inhomogeneous inclusions with e(x) = ju(x) are equivalent to 
empty space and are thus invisible. Pendry, Schurig and Smith |20j went one step fur- 
ther and, similar to what Greenleaf, Lassas and Uhlmann [21] had done in the context 
of the conductivity equations, found that by using a transformation which maps a point 
to a sphere one could create a cloak having e(x) = /x(x) which completely shields an 
object, making it invisible to a fixed frequency of incident radiation. Rigorous mathemat- 
ical justifications of this cloaking, and its generalizations, have been given [22l |23"1 124"] . 
and various approximate cloaks have been proposed and some have been physically con- 
structed [25l [26j |27l EU [29l [30] . This and other types of cloaking are surveyed in reviews 
[3T1 [32] . The growing field of transformation optics uses coordinate transformations to 
achieve unusual effects, such as rotators [33], concentrators |34j, wormholes [35], novel 
superlenses [361 EH EB1 [39] but requires one to tailor materials with prescribed tensors s 
and \x. In the geometric optics limit one can use transformations of the refractive index 
to achieve cloaking [4*01 |4"T] and novel Eaton lenses [12] . 

Extreme values of e and /x, near zero or infinity, also have striking applications, 
leading to new types of circuits [121 HH US] and tunneling of energy through narrow 
channels [4"6] . 

Despite these tantalizing results, many of which are summarized in the book of Cai 
and Shalaev [47], to my knowledge no one has either experimentally, numerically or the- 
oretically shown that any given pair of prescribed symmetric real tensors e and /x, and in 
particular given tensor pairs with e = n, including isotropic materials with e — fi — — 1, 
can be almost realized at a fixed frequency. (Isotropic materials with negative real per- 
mittivity and permeability have been realized [H] , but it is unclear if they can be realized 
with e = fi = — 1.) This paper accomplishes that goal from a theoretical perspective. It 
shows there are no hidden constraints restricting the possible real pairs (e, n) that can 
exist at one given frequency. The construction uses hierarchical microstructures with 
structure on many widely separated length scales, and component metamaterials with 
arbitrarily low loss. It seems likely that the same effective properties could be achieved 
with more realistic designs, but it is doubtful that their properties could be calculated 
analytically. 
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The starting materials and Tartar's formulae for the 
effective tensors of laminates 



In this paper we assume an idealized world where the continuum, and Maxwell's equa- 
tions, extend down to arbitrarily small length scales and we assume the existence of two 
metamaterials. Material A has arbitrarily small microstructure and an effective permit- 
tivity tensor with a bounded real part (s A )' which is diagonal, with exactly one negative 
and two positive diagonal elements and an effective permeability tensor with a real part 
which is not necessarily diagonal but bounded and strictly positive definite, i.e. 
there exist constants /3 A > &a > such that 

P A I > (fJL A )' > a A l. (I ) 

The imaginary parts (s A )" and of these two tensors are assumed to be vanishing 

small, but non-zero, i.e. there exist constants 77 > v > such that 

77I > (e A )" > ul, rjl > (//)" > vl, (2) 

where the loss parameters 77 and v are vanishingly small. [Strictly speaking we should 
talk about a sequence of materials A, parameterized by rj and consider what happens 
in the limit as r\ — > 0, when (e A )' and (fJ. A )' do not depend on 77, while v and the scale 
and material constants of the microstructure within the metamaterial, do depend on 77.] 
Material B has arbitrarily small microstructure and an an effective permittivity tensor 
with a real part {s B )' which is strictly positive definite, i.e. satisfying the bounds 

P B I > (e B )' > a B l, (3) 

for some choice of (3b > a B > 0, and with an effective permeability with a bounded 
real part (/J. B )' which is diagonal, with exactly one negative and two positive diagonal 
elements. The imaginary parts {s B )" and (/x B )" of these two tensors are assumed to be 
vanishing small but non-zero, i.e. to satisfy the bounds 

77I > (e B )" > vl, 77I > (ii B )" > ul, (4) 

where, without loss of generality 77 and u are the same parameters as appear in (J2]). 
[Again strictly speaking we should consider a sequence of materials B, parameterized by 
77 and consider what happens in the limit as 77 — > 0.] 

Material A could be a metamaterial comprised of a cubic lattice of well separated 
cubes, where each cube has a microstructure of highly conducting rods aligned in the x\ 
direction. A rigorous mathematical proof that such a material has the desired effective 
properties, with a more precise description of the needed microgeometry, has been given 
by Bouchitte and Bourel jl9]. Based on the results of Pendry et.al. [50] and Bouchitte 
and Felbacq [51] one might think that an periodic array of highly conducting thin rods 
aligned in the x\ direction would suffice for material A. However, Bouchitte and Felbacq 
[52] show that in a finite sample of such a material the associated effective equations can 
have a non-local behavior. 

Material B could be a metamaterial comprised of a periodic lattice of highly con- 
ducting split ring resonators with axes aligned in the x\ direction with one rod per unit 
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cell. The split rings behave like polarizable magnetic dipoles and if one is just above 
resonance these can have a negative permeability in the X\ direction. This means of cre- 
ating artificial magnetism was discussed by Schelkunoff and Friis [53J , and although their 
formula showed a negative permeability in the X\ direction, they did not draw attention 
to this fact. The microstructure was rediscovered by Pendry et.al. [M], who realized the 
significance of the negative permeability. A rigorous mathematical proof that such a ma- 
terial has the desired effective properties, has been given by Bouchitte and Schweizer |55j, 
following earlier work of [561 EU EHJ [59] [60] (see also the introduction of |61j). Bouchitte 
and Schweizer also prove that one can obtain isotropic effective magnetic permeability 
tensors with a negative permeability. 

We explore what effective permittivity and permeability tensors can be generated from 
composites of these two metamaterials A and B and all rotations of them within classical 
homogenization in the limit as rj tends to zero. The first step will be to recover the result 
of Bouchitte and Bourel [IH] that a material with any real symmetric permittivity tensor 

can be approximately achieved. Then the same argument applies to show that a 
material with any real symmetric permeability tensor /z* can be approximately achieved. 
Finally by combining these two results we show that any pair of real symmetric tensors 
(e*,/zj can be approximately achieved. 

The composites we consider are multiple rank laminates. This class of composite was 
first introduced by Maxwell [62] and has been extensively studied in the homogenization 
literature: see, for example, [63] and references therein. They are obtained by laminating 
the starting (rank 0) materials together on an extremely small length scale to obtain 
rank 1 laminates and then laminating these rank 1 laminates with other rank 1 or rank 
laminates on a much larger, but still extremely small length scale, using a possibly 
different direction of lamination, to obtain rank 2 laminates and so forth. The advantage 
of considering this class of composites is that their effective tensors can be explicitly 
calculated, in the limit in which the ratio of successive length scales approaches infinity. 
The bounds ([2]) and (jlj) ensure that one can use reiterated homogenization theory: to 
calculate the effective tensor in the limit of widely separated scales of a multiple rank 
laminate which is a simple laminate of two (possibly multiple) rank laminates C and D 
one can replace C and D by homogeneous materials with material tensors the same as 
the effective tensors of C and D. 

The use of classical homogenization, except inside the metamaterials A and B, should 
be valid provided, when we replace materials A and B by homogeneous materials with 
material tensors the same as the effective tensors of A and B, the wavelengths and atten- 
uation lengths within each substructure should be much larger than the microstructure 
at that level. For any fixed rj > this should be ensured by ([2]) and (j4|), in the limit 
in which the overall microstructure tends to zero, and the ratio between scales tends to 
infinity. 

If a set of materials are laminated in direction x\ then the formulae obtained by 
Tartar [M] (using an idea that goes back to Backus [65J) for the effective permittivity 
and permeability are 



where s(xi) and /J-{xi) are the local permittivity and permeability tensors (that could 
themselves be effective tensors) which only depend on x\ and the angular brackets denote 
volume averages and for any symmetric matrix C with elements Cjj the matrix C is 




(5) 
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symmetric and has elements 

en = -1/cn, c lk = c lk /c u , c ke = c kl - c kl c u /cn, (6) 

for all k ^ 1, £ =^ 1. Conversely if a symmetric matrix C is given, then C is symmetric 
with elements 

Cll = — 1/Cii, Cife = — Cifc/Cn, C k £ — C k ^ — C k iCx£/Cxi, (7) 

for aU jfc ^ 1, ^ 1. 

In the case where e is diagonal, (jSJ) and (EJ) reduce to the familiar harmonic and 
arithmetic averages 

e* kk = {e kk ) k^l, (8) 

with similar formulae applying when /x is diagonal. 

In the limit as the imaginary part of e approaches zero then (jSJ) implies that the 
imaginary part of £* will also tend to zero unless one is at resonance where the real part 
of (I/en) or the real part of (l//xn) vanishes. We will always avoid this happening. Then 
the limiting tensors e* and /x* can be obtained by replacing in (EJ) £ and n with their 
real parts, i.e. the vanishingly small imaginary parts do not effect the effective tensors 
except at resonance. From now on we will drop the primes and treat e and n as real, 
while remembering that they do have vanishingly small imaginary parts to ensure that 
we can apply the rules of classical reiterated homogenization. 

Realizing a material with a desired tensor e* using 
lamination 

Now let us explore what effective permittivity tensors e* can be obtained by hierarchically 
laminating material A with itself (and its rotations by 90° about the axes) in directions 
parallel to the axes, calculating the effective tensors at each stage using the harmonic and 
arithmetic averages flS}. Let us assume the axes have been chosen so £ A =Diag[— a, b, c], 
with a, b and c all positive. The following remark is helpful: 

Remark 1 If e=Diag[— a, b, c], with a and b both positive is realizable then so is 
£*=Diag[— a/S, bS, c] for any finite 5^0. 

To prove this remark, let us laminate Diag[— a, b, c] with the rotated material Diag[6, —a, c] 
in direction x\ in volume fractions f\ and /*2 = 1 — j\. The resulting material has ef- 
fective tensor e*=Diag[— a/S, bS, c] with S = f\ — f 2 (a/b) taking all values inside the 
interval between 1 and —a/b as fx ranges between and 1, excepting 5 = where one 
of the eigenvalues of s* becomes infinite and one is at resonance. Alternatively, let us 
laminate these two materials in direction The resulting material then has effective 
tensor e*=Diag[— a/ S, bS, c] with S = [fx — }'2{b/ a)] -1 taking all values outside the interval 
between 1 and —a/b as fx ranges between and 1. Using one of the two constructions 
any finite value of S ^ is possible. 

Starting from e A =Diag[— a, b, c], with a, b and c all positive and applying Remark 
1, we obtain a material with effective tensor Diag[a/<5 , —bS ,c], with S > 0. Applying 



5 



remark 1 again, we obtain a material with effective tensor Diag[a/<5o, —bSo/Si, c5\], and by 
rotation (and replacing 8 X with I/S2) a material with effective tensor Diag[a/ S , cj 8 2 , — 65 2 5q] 
Laminating these last two materials together in equal proportions in direction x\ gives a 
material with effective tensor Diag[a/5o, e, g], where 

e = (c/62 - b5o/8 x )/2, g = (c8 x - b5 2 8 )/2. (9) 

Given prescribed values of e > and g ^ we may choose non-zero 

g± ^g 2 - {5 bcg/e) c5 1 -2g 
c bd 

so that OH]) is satisfied provided we choose 5o with 

0<5 < \ge\/(bc), (11) 

to ensure that the roots of (flOl) are real. Since (TTTT) remains valid if we change the sign 
of e we can also realize a material with effective tensor Diag[a/5o, —e,g], and hence by 
remark 1, a material with effective tensor Diag[— a/5o, e, g]. Laminating this material 
in direction x 2 with Diag[a/5o, e, g] we can obtain a material with a prescribed effective 
tensor Diag[/i, e, g] with e > provided 5 is chosen sufficiently small to satisfy (ITT]) 
and the constraint that 5o < a\h\ which ensures h lies between a/5 and —a/5o- So 
any diagonal tensor with non-zero diagonal elements, at least one of which is positive, is 
realizable. 

The only case left to treat is when h, e and g are all negative. Let t be bigger that 
both — e and —g. Then the tensors Diag[/i, e + t,g — t] and Diag[/i, e — + are 
both realizable and by laminating them together in direction x x we see that the tensor 
Diag[/i, e, g] is also realizable as an effective permittivity tensor for any non-zero finite 
choice of h, e and g. 



Realizing any desired pair of tensors 

By rotating the material just obtained we see that any symmetric matrix e A *, except 
possibly those which are singular, is realizable as an effective permittivity tensor. The 
associated effective permeability tensor fj, A * (in the limit the loss goes to zero) must 
satisfy the classical homogenization Weiner bounds (66] 

<M> > V M > (12) 

where /x(x) is the local permeability tensor, which is locally a rotation of n A . Since \x A 
satisfies ([TJ we conclude that /j, a * also satisfies the inequality 

P A 1 > fi A * > a A I. (13) 

Let us call U the family of materials thus obtained. 

Similarly (because we can interchange the roles of e and fi) using material B, and 
its rotations, any symmetric matrix /j, b *, except possibly those which are singular, is 
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realizable as an effective permeability tensor, and the associated effective permittivity 
e B * satisfies the bounds 

PbI > e B * > a B I. (14) 

We let V denote this family of materials. 

Now we are free to take any set of materials in the set W = U U V and laminate 
them together in direction x% to form a larger set of materials LW containing W. By 
an abuse of notation we let W also denote the set of tensor pairs (e*, Jl*) deriving from 
effective tensors (e*, /x*) of materials in W, and we let LW also denote the set of tensor 
pairs (?*, Jl*) deriving from effective tensors (e*, /j,*) of materials in LW. Now fl5]) implies 
the tensor pairs in LW correspond to all convex combinations of the tensor pairs in W, 
except for possibly those tensor pairs for which either e\ x = or Ji*^ = where one is at 
resonance. So aside from a set of measure zero the set of tensor pairs LW is a convex set 
whose convex hull is characterized by the Legendre transform: 

3 3 

/(P, Q) = inf + ( 15 ) 

as a function of the symmetric matrices P and Q. If we can show that f(P, Q) is minus 
infinity for every choice of P and Q, not both zero, then we can conclude that pair of 
real symmetric matrices (e*,fl*) is approximately realizable. 
First observe that 

/(P, Q) = min{/ A (P, Q),/ B (P, Q)}, (16) 

where 

3 3 

/a(P, Q) = inf Yl XX^* + ^S*' 
(£ *)eu i=1 j=1 

3 3 

/b(P,Q)= B inf ^»/:f • 'luf'u- (17) 
(£ *,M *)ev i=1 i= i 

Since almost every matrix e A * is realizable as a permittivity tensor, it follows that almost 
every matrix e A * is realizable: given e A * we realize, or almost realize, the tensor e A * whose 
elements are given according to ([7]). In particular we can realize a tensor e A * — AP, 
for any value of A no matter how large. The associated permeability tensor fi A * satisfies 
(TT5|) . which implies 

(3 A > /if > a A , f3 A > |^*| for all i ^ j, (18) 

where the latter condition follows from the positivity of the determinant of the associated 
2x2 submatrix of [i A *. It follows that the elements of fl A *, given according to fl6]), satisfy 
the bounds 

lA*n I < 1 / 0£ a, \Jitk \ < P>aI^a, I Mm I < Pa + Pa/ola, (19) 
irrespective of the choice of A. Consequently, for any fixed choice of P ^ and Q, 

i=l j=l 
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approaches — oc as A — > oo, while 

3 3 

i=i j=i 

remains bounded. It follows that /a(P, Q) is minus infinity unless P = 0. Similarly 
/b(P, Q) is zero unless Q = and thus /(P, Q) is minus infinity unless both P and Q 
are zero. We conclude that any given pair of real symmetric matrices (e*, fj,*) is approxi- 
mately realizable as the (effective permittivity, effective permeability) of a metamaterial. 
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